A Mathematica interface to NNPDFs 



Nathan P. HartlancP, Emanuele R. Nocera b 

"Tait Institute, School of Physics and Astronomy, University of Edinburgh, EH9 3JZ, UK. 
1 ' Dipartimento di Fisica, Universita di Milano and INFN, Sezione di Milano, Via Celoria 16 1-20133 Milano, Italy 



Abstract 

We present a Mathematica interface for handling the parton distribution functions of the NNDPF Collaboration, 
available from the NNPDF hepforge website http : //nnpdf . hepf orge . org/| As a case study we briefly summarise 
the first PDF set which includes all relevant LHC data, NNPDF2.3, and demonstrate the use of our new Mathematica 
interface. 
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We have developed a Mathematica package which 
provides direct interactive access to any NNPDF par- 
ton set. This package reads the NNPDF LHgrid files 
as available from the LHAPDF library and performs 
their interpolation in (x,Q 2 ) space by means of proper 
Mathematica built-in functions. The user can there- 
fore make use of PDF ensembles within a Mathemat- 
ica notebook in order to compute observables, evaluate 
PDF central values and variances and plot PDFs instan- 
taneously. In this note we shall present a brief motiva- 
tion and summary of the key features present in the new 
NNPDF2.3 0] PDF set, and demonstrate the NNPDF 
Mathematica package by examining this new determi- 
nation. 

In recent years the volume of collider data of applica- 
tion in the accurate determination of parton distribution 
functions has significantly increased. Previous PDF de- 
terminations based solely upon high energy collider data 
collected by Tevatron and HERA collaborations were 
poorly constrained in comparison to global fits includ- 
ing data from fixed target deep inelastic scattering and 
Drell-Yan experiments. These collider only determina- 
tions are of particular interest in that they are free from 
the potential issues arising from nuclear corrections and 
higher twist effects present in the low energy data of 
global fits. Studying the impact of LHC data upon PDF 
determinations is therefore crucial not only for further 
constraint upon global determinations, but also for their 
potential utility in providing competitive collider only 
fits. 

The NNPDF2.2 parton set [2 1 was the first determina- 
tion to investigate the constraining power of LHC mea- 
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surements upon parton distributions, specifically the 
W lepton charge asymmetry measurements of the AT- 
LAS and CMS collaborations in addition to the 
DO collaboration measurements |5] [6). This data was 
included by taking advantage of the Monte Carlo uncer- 
tainty estimation in the NNPDF methodology. NNPDF 
sets provide a probability distribution in the space of 
PDFs allowing for the use of a Bayesian reweighting 
technique [7, 8 | and the rapid inclusion of new datasets 
into a prior set. However, this method is limited to the 
inclusion of a relatively small dataset. 

In order to include an updated and more compre- 
hensive LHC dataset into an NNPDF fit, an extremely 
fast method of collider observable computation has been 
developed. The FastKernel method for hadronic ob- 
servables [ 1 ] combines the Monte Carlo weight grids 
from fast NLO software packages such as APPLgrid Q 
and FastNLO JTOJl, with FastKernel PDF evolution ta- 
bles BllL The resulting FastKernel table may then be 
combined with the initial scale parton distributions to 
produce the required observable. For a given hadronic 
observable cr the operation required to compute the the- 
oretical prediction for such an observable is simply, 

-Wpdf N x 

where 2 is the final FastKernel table and the N are the 
initial scale, DGLAP evolution basis parton distribu- 
tions. The indices i, j run over the N x points in the fitting 
scale x-grid, and a, f3 run over the reduced flavour basis 
of A^pdf light PDFs at this initial scale. This operation 
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is computationally very simple and allows for a much 
accelerated fit. 

With this method we are able to add a relatively large 
LHC dataset to the existing data in the NNPDF 2. 1 fam- 
ily to produce a set of fits denoted NNPDF 2.3. The 
NNPDF2.3 LHC dataset includes ATLAS collaboration 
35 pb~' measurements of the inclusive jet cross sec- 
tion lfl2l and electroweak vector boson rapidity distribu- 
tions H131 . The CMS collaboration's 840 pb~' measure- 
ment of the W electron asymmetry [ 14] and the LHCb 
W boson rapidity distribution data 1TT31 are included 
also. The expanded kinematic reach of the NNPDF2.3 
dataset is described in Fig.[T] In addition to the inclusion 
of the LHC data, the NNPDF2.3 family of fits features 
a number of methodological improvements. 



NNPDF2.3 dataset 




Figure 1: Kinematic range of the NNPDF2.3 dataset. 

In order to assess the impact of the included LHC 
data, we compare the NLO NNPDF2.3 [ 1 1 and NLO 
NNPDF2.1 02) singlet sector PDFs at Q 2 Q = 2 GeV 2 
(see Fig [2]). We have used our new Mathematica inter- 
face to draw these plots. From Fig. [2j it is clear that 
PDFs from the two sets differ by less, and usually much 
less, than one sigma. 

The main advantages of our interface come from the 
possibility to combine the use of any NNPDF set with 
all the Mathematica features. In particular, this al- 
lows to perform a variety of manipulations on PDFs 
in a straightforward and interactive way, as we briefly 
demonstrate in the following discussion. 

• Compute PDF central value and variance We 

have defined proper functions to keep the compu- 
tation of PDF central value and variance very easy. 
These built-in functions only need x, Q 2 and PDF 
flavour as input. The user can also specify the con- 
fidence level to which central value and variance 
should be computed. 
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Figure 2: NLO NNPDF2.3 (green) singlet sector PDFs at Q\ = 
2 GeV 2 , compared to their NNPDF2.1 (red) counterparts. All 
error bands shown correspond to a one sigma interval. These 
plots are drawn with the new NNPDF Mathematica interface. 
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Figure 3: Simultaneous dependence from both x and Q 2 for the 
one sigma error of the NLO NNPDF2.3 Singlet PDF. This plot 
is drawn with the new NNPDF Mathematica interface. 



• Make PDF plots Mathematica enables a wide 
range of plotting options. As a few examples, we 
show the 3D plot and the contour plot of the Sin- 
glet PDF from the NLO NNPDF2.3 parton set (see 
Figs. [3]-|4]respectively). 

• Perform computations involving PDFs PDF ma- 
nipulation can be carried out straighforwardly 
since we have defined functions which handle ei- 
ther single replicas or the whole Monte Carlo en- 
semble. The user can then easily perform any com- 
putation which involves PDFs. For example, we 
show in Fig.|5]a snapshot of a typical Mathematica 
notebook in which we use our interface to check 
the momentum and valence sum rules from the 
NLO NNPDF2.3 parton set. 

The new NNPDF Mathematica package can be 
downloaded from the NNPDF web site, 



http : / / nnpdf . hepf orge .org/ 



together with a sample notebook containing a step by 
step explanation of the NNPDF usage within Mathemat- 
ica as well as a variety of examples. 
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Figure 4: Contour plot of the NLO NNPDF2.3 Singlet PDF in 
the (x,Q 2 ) plane. This plot is drawn with the new NNPDF 
Mathematica interface. 
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Figure 5: A snapshot of the Mathematica notebook written with our interface for checking the momentum and valence sum rules 
from the NLO NNPDF2.3 parton set. 
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